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Abstract. We introduce a geometric evolution equation of hyperbolic type, 
which governs the evolution of a hypersurface moving in the direction of its 
mean curvature vector. The flow stems from a geometrically natural action 
containing kinetic and internal energy terms. As the mean curvature of the 
hypersurface is the main driving factor, we refer to this model as the hyperbolic 
mean curvature Bow (HMCF). The case that the initial velocity field is normal 
to the hypersurface is of particular interest: this property is preserved during 
the evolution and gives rise to a comparatively simpler evolution equation. 
We also consider the case where the manifold can be viewed as a graph over 
a fixed manifold. Our main results are as follows. First, we derive several 
balance laws satisfied by the hypersurface during the evolution. Second, we 
establish that the initial-value problem is locally well-posed in Sobolev spaces; 
this is achieved by exhibiting a convexity property satisfied by the energy 
density which is naturally associated with the flow. Third, we provide some 
criteria ensuring that the flow will blow-up in finite time. Fourth, in the case 
of graphs, we introduce a concept of weak solutions suitably restricted by 
an entropy inequality, and we prove that a classical solution is unique in the 
larger class of entropy solutions. In the special case of one-dimensional graphs, 
a global-in-time existence result is established. 



1. Introduction 

Our aim in this paper is to introduce and study a geometric evolution equation 
of hyperbolic type which describes the flow 

F : [0, T) x M — ► M" +1 , T > 

of an immersed n-dimensional hypersurface M in the Euclidean space. We derive 
this evolution equation from a geometrically natural action functional based on the 
local energy density 



dt 



involving the kinetic energy of the hypersurface and the internal energy associated 
with its volume. The equation under consideration models the nonlinear motion of 
an elastic membrane, driven by its surface tension only. Our model is purely geo- 
metric and requires no constitutive equation on the membrane material (contrary 
to what is required in the theory of nonlinear elastic bodies or shells) . As the mean 
curvature of the hypersurface is the main driving factor, we refer to this model as 
the hyperbolic mean curvature flow (HMCF); see Proposition ^. 2l below. Stationary 
solutions of this flow will be minimal hypersurfaces with vanishing kinetic energy. 
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The flow equation takes a simpler form in the case that the initial velocity is 
normal to the hypersurface, i.e. if its tangential part vanishes: 



dt 



= 0. (1.1) 

t=o 



Namely from the momentum conservation law satisfied by a general flow, we will 
deduce that tangential components of the velocity vector vanish for all times if 
they vanish initially. Hence, under this assumption, the (normalized version) of the 
proposed HMCF equation reads 

d2F rr r, 

t (HMCF') 
= 0. 



dF 
dt 



\t=o 



where the scalar H is the mean curvature of the hypersurface and the vector v 
denotes its unit normal (chosen to be inward pointing when M is compact without 
boundary). In fact, the assumption (jl.ip is geometrically motivated in the sense 
that tangential variations do not alter the shape of the hypersurface and merely 
correspond to reparametrizations by a suitably chosen family of (time-dependent) 
diffcomorphisms. Since, geometrically, only (HMCF') is of interest, we will mainly 
study this flow, which we refer to as the normal mean curvature flow equation. 

The main results established in the present paper are as follows. After introduc- 
ing the proposed flow in Sections [2] and [3j we derive in Section[4]several conservation 
laws or balance laws satisfied by the hyperbolic mean curvature flow. Then, in Sec- 
tion [5l we begin our investigation of the properties satisfied by general solutions to 
the hyperbolic flow by restricting attention to the important case that the hyper- 
surface is represented as an entire graph over 1™ : we prove the local well-posedness 
of the flow equation, and introduce a concept of weak solutions suitably restricted 
by an entropy inequality; we also prove the uniqueness of a classical solution within 
the class of weak solutions, and for one-dimensional graphs we establish the global- 
in-time existence of weak solutions with bounded variation. The convexity of the 
measure e d/i with respect to certain well-chosen variables is an essential observation 
for these results. Then, for the rest of the paper we turn to the equation (HMCF') 
for normal flows and, in Section [6l we prove that the equations under consideration 
can be recast in the form of a first-order nonlinear hyperbolic system, and we ob- 
tain a local- in-time existence result for the evolution of general compact manifolds. 
Next, in Section [3 we provide some criteria ensuring that the flow will blow-up 
in finite time, due to the formation of geometric singularities or shock waves. For 
general material on flows by mean curvature we may refer to [1] [3] , and on nonlinear 
wave equations to [2 [5j [6] . 



2. Structure equations for general flows 

Let F : M — > R n+1 be a smooth immersion of an orientable smooth manifold 
M of dimension n into and let v be the unit normal vector defined along 

the hypersurface and chosen to be inward pointing when the manifold is compact 
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without boundary. In local coordinates {x l )i = i n , we have 
and 

F x A . . . A F n 



^a.^af^ 

The induced metric g = g^ dx % ®dx^ and the second fundamental form h — dx l ® 
dx J of the hypersurface are 

9ij = ( F i > F j)i h H =-( F i^3 v )i 
respectively. Here, V denotes the Levi-Civita connection associated with g. Through- 
out, we use Einstein's summation convention on repeated indices and, for simplicity, 
we keep the same notation (•, •) for both the standard inner product on R™ +1 and 
the induced inner product on M. Latin indices are raised with the inverse (</ y ') of 
the metric (gij) so, for instance, 

h{ : h ik g kj . 

We denote by Rijki the components of the Riemann curvature tensor of the 
hypersurface in local coordinates. We denote the induced volume form on M by 
d\x and the (scalar) mean curvature by H := g^hij. We will use also the following 
convention: we identify the gradient Vp of a function p on M with its image 
dF(Vp) = V'pF,. 

The following basic properties of these tensor fields are easily checked from their 
definitions. 

Lemma 2.1. The Gauss- Weingarten-Codazzi equations of the hypersurface M 
read: 

WiFj = hijis, (2.1) 

Vns = -h{Fj, (2.2) 

Viftjfc = Vjh ik , (2.3) 

Rijki = hikhji - huhjk- (2.4) 

In the present paper, we are interested in a How of hypersur faces, that is, a 
smooth family of immersions 

F : [0,T) x M -> R n+1 , 

so that all of the tensor fields defined above also depend on the time variable t. We 
can then define on M some additional (time-dependent) functions a, a and 1-form 
fields S = Sidx\A = Adx 1 by 

d -F,J\, Si:=/^FFi 



dt ' / ' 1 ' \ dt 

d* \ Id 2 



We refer to a and S as the normal and tangential velocity components, respectively, 
and to a and A as the normal and tangential acceleration components, respectively. 
We have the decomposition 

^-F = au + S i F i , f- F = a v + A*F l . 
dt dt 2 
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To express the structure equations satisfied by a general flow, it is convenient to 
introduce the local energy density of the hypersurface 



1 

G:= 2 



dt 



d 
dt 



(2.5) 



where we have used VF = g lJ (VjF, VjF) = g v gij — n. A simple computation 
based on Lemma 12.11 then yields the following expressions for the components of 
the velocity and acceleration fields. 

Lemma 2.2. Every flow F satisfies the following structure equations: 

V i a=(^F i ,v\-hP i S j , (2.6) 



dt 
d 

di~ 



ViS^i-rFiiFA+o-hi^ (2.7) 



a = j t a+ (V(T,S) +h(S,S), (2.8) 

A = 4- S - de, (2.9) 
dt 

where de = J \7jedx J denotes the exterior differential of the function e, dual to the 
gradient Ve = V l eFi. 

The following notion will be of special interest in this paper. 

Definition 2.3. A flow F : [0, T) x M R" +1 is called a normal Row if and only 
if its tangential velocity vanishes identically, that is, S{t) = for all t € [0, T). 
Proposition 2.4. 

1. A flow F is normal if and only if its tangential velocity and tangential accelera- 
tion satisfy 

5(0) = at the initial time, 
A = — de at all times. 

2. Given a general flow F : [0, T) x M — > M™ + , there always exists a smooth 
family of time- dependent diffeomorphisms : M — > M such that the modified 
flow given by 

F : [0,T) x M -> R" +1 , F(t,x) := F(t,^ t (x)) 

is a normal flow. In particular, the hypersurfaces 

M t :=F{t,M), M t :— F(t, M) 

coincide for each t £ [0,T). 



Proof. The first claim follows immediately from equation (|2.9p . To derive the sec- 
ond claim we consider a general flow 

— F(t, x) = a(t, x) u(t, x) + S(t, x). 
dt 



THE HYPERBOLIC MEAN CURVATURE FLOW 



5 



Since S(t,x) is tangential to Mt, we can introduce the solution '■ M — » M be 
the following ordinary differential equation (ODE) 

^* t (a:)=-S(t,* t (a:)). 
Then, by setting S(t,x) =: S l (t,x) Fi(t,x) we see that the map 

F(t,x) :=F(t,* t (x)) 
satisfies the evolution equation 

j t F(t, x) = a(t, * t (i)) u(t, M*)) + S(t, M*)) + F i(t, Mx))j t *i(x) 
= a(t,V t (x))v(t,y t (x)) 
=: a(t, x) u(t, x). 

□ 

The proposition above shows that a general flow F and its normalized version F 
can be identified geometrically. Therefore, without genuine restriction, our analysis 
will often be focused on normal flows, which have the general form 

4~F = av -Ve. (2.10) 
at 1 

It should be observed that, at this stage, the normal acceleration a has not been 
defined yet. Our results will show that by prescribing this scalar field (in the 
forthcoming section) the evolution of the hypersurface is uniquely determined. 

To conclude this section, in view of the computations done in Huisken [3] for 
the standard mean curvature flow, we obtain the following first-order evolution 
equations for the induced metric, volume form, second fundamental form, and mean 
curvature of the hypersurface. 

Lemma 2.5. The evolution of the tensor fields gij^dn^v^hij,!! associated with a 
general flow F : [0, T) x M — -> M™ +1 is determined by the equations 

— gij = -2ahij + V l S j + VjSt, (2.11) 

—du=((FS-o-H)du,, (2.12) 
dt 

^-v = -{V i a + h ik S k )F u (2.13) 
dt 

£hij = ViVj-ff - trhfhkj + h^VjSk + h%ViS k + V k h i:j S k , (2.14) 

^-H = Aa + a\h\ 2 + S l ViH, (2.15) 
dt 

where <v S — W l Si — WiS 1 denotes the divergence of a vector field, and \h\ 2 := h l i hij 
denotes the (squared) norm of a 2-tensor field. 
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In particular, if the flow is normal we take S 
and obtain 



and A 



d 

dt 9ij ' 



— d/i = —crHdu, 
dt 

d 

-u=-V l aF l , 
dt 

d 

dt 
d 

dt 



hij = ViVjfT — crh i hkj, 
H = Aa + a\h\ 2 . 
3. The hyperbolic mean curvature flow 



de in Lemma 12.51 

(2.16) 
(2.17) 
(2.18) 
(2.19) 
(2.20) 



We are now in a position to introduce the evolution equation that we propose in 
this paper. The flow is going to be defined from an Hamiltonian principle based on 
a geometrically natural action, consisting of a kinetic term and an internal energy 
term, which is defined geometrically as the local volume density of the hypersurface. 

More precisely, let F : [0,T] x M — > M' i+1 be a smooth family of immersions of 
an orientable manifold M of dimension n into R' i+1 . Define the total kinetic energy 
at the time t by 



K(t) 



1 



M 



dt 



-F 



d[i t 



and, after integrating over the time interval [0, T], consider the action 



dfitdt. 



Define the total internal energy of the hypersurface at the time t by 







f I - 




'o Jm 2 


dt 



(3.1) 



V(t) := 



1 



|VF| 2 ^ t - 



A I 



M 



solely determined by the induced volume form, and consider the corresponding 
action 



Jv(F) := J / / dfitdt. 
z Jo Jm 



(3.2) 



According to the Hamiltonian principle, we impose that the evolution of the 
hypersurface is stationary for the action Jy — Jk, that is, 

d 



ds 



(J v ~ J K ){F + s$) 



s=0 







(3.3) 



for all$ G Cg°([0,T]xAf,] 
oat any order). 



(compactly supported maps that are differentiable 



Remark 3.1. Obviously, if the volume of the manifold is infinite, the functionals 
Jk{F) and Jv{F) above are only formally defined. This difficulty can be easily 
overcome by restricting attention to any compact subset of M. However, since the 
stationarity condition (|3.3[) implied by the Hamiltonian principle itself is formulated 
in terms of compactly supported variations, this is unnecessary. 



We now show: 
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Proposition 3.2 (Hyperbolic mean curvature flow equation). The stationary so- 
lutions of the action functional Jy — Jk satisfy the equation of motion 

^F=au + A k F k , 

' a:=(e-\S\ 2 )H-atfS, (HMCF) 

A k :={(jH - SS)S k - V k e. 

For instance, if the initial velocity is normal, i.e. if 6*14=0 = 0, then it will follow 
from the conservation of momentum that Su = for all t, so that the HMCF 
equation reduces to the much simpler equation 

d 2 

—F = eHv-Ve, 

T (HMCF') 

—F ) = 0. 

dt /|t=o 

Comparing with the general flow equation (|2.10|) we see that the minimal action 
principle allows us to identify the normal acceleration, as a linear function in the 
mean curvature H: 

a = eH. 

In consequence, for normal flows the evolution of the normal component of the 
velocity is proportional to the mean curvature 

Later in this paper, we will prove that the equation (HMCF') is hyperbolic. We 
will not treat here the general system (HMCF) except in the next section where we 
will derive some (a-priori) conservation laws for the general flow. The subsequent 
sections are entirely devoted to the normal equation (HMCF'), since this (as already 
noted in the introduction) is sufficient from the geometric point of view. 

Observe also that Ve depends on second-order derivatives of F, namely on mixed 
derivatives in space and time. In agreement with the standard mean curvature flow 
(which is parabolic) , the acceleration a is defined in terms of the curvature of the 
manifold M, and a, considered as an operator on M, is elliptic. Therefore, this 
justifies to refer to the proposed flow as the hyperbolic mean curvature Row. The 
following sections will show, both, some similarities and some marked differences 
between the parabolic and hyperbolic versions of the mean curvature flow. 



Proof. A simple computation yields 

d_ 

ds 



d J v (F + s$)\ s=Q = ^ f j ^{F^d^dt 

z JO JM 



[ H{u, $)du. t dt, 

JM 



n 

where, for the second identity, we have integrated by parts and used the contracted 
Gauss formula AF — Hv (a consequence of (|2.Ijl ). 

On the other hand, for Jk we obtain the first variation formula 

±MF + »*)!„„= fj M I*) + (e- ^ 9 >HF„* l} ) 
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in which we now successively integrate by parts each term of the right-hand side. 
For the first term we find 



/q \ dt J dt 



l L(-^ F+iaH - dfs) ^) d ^ 



where we used (|2.12l) . For the second term we obtain 

i-T r r T 



I I (e- ~)g ij (F i ,$ j )di H dt = - [ [ <Ve+ (e- ~)Hp,$)d» t dt. 
Ja Jm 1 Jo Jm 1 



Combining the above identities together, we deduce that 

r-T 



^-(Jv~Jk)(F + s<P)\ s=0 = f f (P,Z>)d» t dt 
as Jo Jm 



with 



d 2 d 
P :=—,F-{aH-d)S)—F + Ve+(e-n)Hv 
at' 1 dt 

=^F + Q(|Sf -o 2 - n)H + ad^S^j v + (V fc e - (aH - d f S)S k ) F k . 

Since <I> S Cq°([0,T] x M, R n+1 ) is arbitrary this completes the derivation of 
(HMCF). □ 

4. Conservation laws and balance laws 

In this section we derive various conservation laws satisfied by solutions of the 
hyperbolic mean curvature flow (HMCF), and we show that, as far as the geometry 
of the hypersurface is concerned, one may work within the class of normal flows 
(HMCF'). Consider a family of immersions F : [0,T) x M -> W l+1 satisfying the 
hyperbolic mean curvature equation (HMCF). In the present section, all tensor 
fields under consideration are assumed to be sufficiently smooth. We use the so- 
called abc method discussed by Shatah and Struwe in [5], and multiply (HMCF) 
by an expression of the general form 

a^-F + b-VF + cF, 
dt 

by choosing the variable coefficients a, b, and c so that higher-order terms in the 
corresponding evolution equations admit a divergence form. 

The following lemma shows that there exists a divergence-type form for (HMCF), 
which will be useful for the derivation of conservation laws (modulo lower-order 
terms). 
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Lemma 4.1 (A general identity). Let Y — Y(t, y) be a time dependent vector field 
on K n+1 . Every solution of (HMCF) satisfies 

J t ((j t F - Y ) d l") = h"((n-em,Y)) + (e-n)g«DY(F i ,F j ) 

+DY {l F -i F ) + {i F - Y ))^ <«) 

where Y t = -^Y and DY is the spatial differential ofY, i.e. 

dY 1 

DY = 5 ai -^dy a ®dy t3 . 
Proof. We multiply (HMCF) by Y and compute 

((e- \S\ 2 )H - arf S\(v,Y) + UaH - <& S)S k - V fc e) (F k , Y) 

i F ' Y )^ s -^ + DY (i F 4 F ) + (i F - Y -)) 

Since -j^F = av + S k F k we obtain 



1 (l—F,Y 

dt \\dt 



- '(e- lS f-^H(,,n-^,Y } + DY(f t F,lF) + (lF,Y t )) «„„ 
= (in - c)H{v, Y) - V*e(F,,F) + By f ±-F, + (i-F, Y,\) dp. t . 



dt dt ) \dt 
From V k F k = AF = 72" we get 

(n - e) JT(i/, y> - V*e<F fc , Y) = V fe ((n - e)<F fc , Y)) - (n - e)g ij DY{F h Fj). 

Inserting this in the above expression for ^ ((^F, Y) d/it) we arrive at gj}. □ 

From Lemma 14.11 we now derive various conservation laws or balance laws of 
interest. The first result below is a consequence of the invariance of (HMCF) under 
isometries of K™ +1 . 

Proposition 4.2 (Local continuity equation). Every solution of (HMCF) satisfies 
|((| F ' F W) = V fc ((n-e)(F fe ,Y))d^, (4.2) 
where Y = is any time independent Killing vector field on K™ +1 . 

Proof. Let Y = be a Killing vector field on K n+1 , that is, Y generates an 

isometry on R" +1 . Since Y is a Killing vector field, DY is skew-symmetric, so that 
the last three terms in (14. ip vanish and we get (|4.2[) . □ 
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Proposition 4.3 (Local momentum equation). Every solution of (HMCF) satisfies 
the balance law 

Jt {(^t F,F ) dfl ) = V ' £ (("- c )( F '^))^+((« + 2)c-n( ? i + l))d A1 . 



Proof. We apply Lemma 14.11 to the vector field Y(y,t) :— y. Then DY(v,w) — 
(v,w) and Y t = 0. Moreover Y(F(x,t),t) = F(x,t), so that (j4~Tj) becomes 

J t ((j t F ' F ) d ^) = (v fc ((n-e)(F fe ,F))+( C -n)^(F J , J F :? ) 

dt F >dt F )) 

= (V((n-e)(F fc ,F)) + (e - n)n + 2e - n) dm, 
which establishes the desired identity. □ 

We next turn to the component S of the velocity vector. 

Proposition 4.4 (Local tangential velocity equation). Every solution of (HMCF) 
satisfies the conservation law 

|(S(X)d/x t )=0, (4.3) 
where X e T(TM) is any time-independent vector field on M. 

Proof. Given a time-independent vector field X = X l -J~ defined on M, from (|2.9|) 
we obtain 

^( SW )= rf e(I) + (^^)), 

= (aH -tflS)(S{X)), 

where the second identity follows by multiplying (HMCF) by X l Fi = dF{X). The 
desired conclusion is now clear in view of (|2.12[) . □ 

The momentum conservation law has the following important consequence. 

Corollary 4.5 (Reduction to normal flows). Within the class of flows F : [0, T) x 
M —> R ,l+1 whose velocity vector is initially normal to the hyper surf ace, i.e. 

F(0, x) = F (x), j t F(0, x) - /(zMO, x). 

where F : M — > K n+1 is a immersion and f : M — ► R a function, the following 
two properties hold: 

(1) the flow F is a solution to (HMCF) if and only if it is a solution of the 
normal flow equation (HMCF'). 

(2) any solution of (HMCF 7 ) satisfies S = for all t £ [0, T); in other words, 
there exists a family of functions a : [0,T) X M — * R such that cr(0,a;) = 
f(x) and 

j t F(t,x) = a(t,x)v(t,x). (4.4) 
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Proof. For each compactly supported tangent vector field X 6 To(TM), by defining 
the total tangential momentum in the direction X as 

p(t,X) := f S(X)dfH, (4.5) 

it is clear that 

S\ t =0 if and only if p(t, X)=0, X 6 T Q (TM). 
However, the identity f|4.3|) implies 

P (t,x)=KO,x) = o, ie[o,T], 

and 

s\ =o => S\=o, te[o,T}. 

Then, we obtain (HMCF') by inserting S = into (HMCF). □ 

We continue our derivation of conservation laws satisfied by the hyperbolic mean 
curvature flow. 

Proposition 4.6 (Local energy identity). Every solution of (HMCF) satisfies the 
conservation law 

j (edfi) = d T ((n - e)S)dfi. (4.6) 

In particular, if the initial velocity is normal along the hypersurface, then e d/i is 
conserved along the flow, 

|(ed/i)=0 (4.7) 
and ed/i can be seen as a fixed volume form on M . 
Proof. We multiply (HMCF) by dF/dt and obtain 

e = l^-F, ((e - \S\ 2 )H ~ ad) S) v + ((aH - d^S)S k ~ V fc e) F k 



= a((e - \S\ 2 )H - ad^S) - V k eS k + (aH - S)\S\ 
= aHc+ (n - 2e)d^S - V k eS k 



2 



= (aH -(tfS)e + <tf((n-e)S). 
Then, the desired identity again follows from (|2.f 2|) . The second statement follows 



from the reduction principle in Corollary 14.51 □ 

From the above proposition the following global result follows. 

Corollary 4.7 (Global conservation laws). If F : [0, T) x M — > W 1+1 be a solution 

of (HMCF') and M is a compact manifold without boundary, then the total energy 
E(t) defined by 

E M (t) := f edfr = \ f dfi t + ^ f dm 

Jm 1 jm at 1 jm 

is conserved: 

E M (t)=E M (0), te[0,T). 

We also state here a compatibility between derivatives of F in time and in space 
which follows immediately from the fact that the derivatives d/dt and V commute. 
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Proposition 4.8 (Compatibility relation between time and space derivatives). 
Every flow F : [0, T) x M — > R" +1 satisfies the conservation law 

±F-V(c-v + &F j )=0. 

Remark 4.9. There are other globally conserved quantities, which are conserved 
for topological reasons and not because of the special nature of our flow. We give 
two examples. 

1. If Y is a divergence free vector field in R™ +1 and M is closed, then by Stokes 
theorem 

J { V ,Y)dfH = 0. 

M 

2. The quantity 

M 

with p e arbitrary (such that F(t, x) ^ p for all (t, x) G [0, T) X M) is the 

degree of the "winding map " and hence an invariant. 

5. Evolution of entire graphs 

In the rest of this article we restrict our attention to the normal hyperbolic mean 
curvature flow (HMCF'). 

Local well-posedness result for the normal flow of graphs. It is convenient 
to begin our investigation with the case of graphs. That is, in this section we 
discuss the case, where F : [0, T) x M -> M™ +1 satisfies (HMCF') and such that 
each M t := F(t, M) is an entire graph over a flat subspace Z 1 - C M™ +1 , where Z 1 - 
denotes the orthogonal complement of a unit vector Z £ R n+1 . 
Without loss of generality, we assume that F is given by 

F(t,x) - (x{t),u(t,x(t))') 

for a time-dependent family of height functions u : [0,T)xM^l and a family of 
diffeomorphisms x(t) = ^.x 1 (i), . . . , x n (t)j defined on the hyperplane 

M = Z 1 - = {(x\...,x n+1 ) e R n+1 : x n+l =0}. 
In such a situation, a solution of (HMCF') is completely determined by the time- 
dependent function u. From now on, we will use the notation utt '■= %pr and 

Ut 3 ■— dtdxo ■ 

Theorem 5.1 (Local well-posedness for the normal flow of graphs). 
1. The hyperbolic mean curvature flow for graphs over a flat hypersurface Z 1 - takes 
the form of the following second-order hyperbolic equation: 

- u tt + (e{a) g l i + a 2 (g lJ - ^)) Uij + 2 ^ 6 ij Ui u tj = 0, (5.1) 

where 

w := + \Du\ 2 = \Jl + S^UiUj, a — u t /w 

and 

g 13 = 5 ij -w- 2 5 lk 8 3l u kUl . 
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w 1 (5.2) 



2. The HMCF equation for graphs can be recast in the form of a nonlinear hyperbolic 
system of n + 1 equations in the unknowns a and b = Du = (ui)i=i,..., n 

da d_{ , h 

dt dxA e[(T) w 

which, moreover, has a conservative form and is endowed with the mathematical 
entropy function 

E(a, b) = ^(°- 2 + n) y/l + S^bibj. 
Moreover, this function is strictly convex provided 

\Du\ 2 = 5 ij bibj < ~. 

3. Furthermore, the graph equation (|5.1[) is locally well-posed in the following sense: 
given data (a, Du) prescribed at the initial time t — and belonging to the Sobolev 
space H S (M" ) (that is, whose all their s- order derivatives are squared integrable) 
for some s > 1 + n/2, there exists a classical solution 

u: [0,T) xF->l m L°°([0,T),H s+1 (W 1 )) D Lip([0, T),H s (R n )), 

defined on a maximal time interval [0, T). 

Proof. Step 1: Derivation of the HMCF' graph equation. We will use that 

d ( d d \ ( d d d A , r , 

—F = —x, —u = —x, —u + Ui—x = ov. (5.3) 
dt \dt dt J \dt ' dt dt J y ' 

On the other hand, since 

d d 

ox % ou 

we obtain 

u = - (-S ij u — + — 
w \ 1 dxi du 



with w := v/l + \Du\ 2 = yl-j- S^UiUj. Inserting this into (|5.3p gives the two 
equations 

dt w 



-y = --S ij Uj , (5.4) 



d d ■ a 

-^-u + Ui—x = —. (5.5) 
at dt w 



Inserting (|5.4p into (|5.5[) yields 



u = aw. (5-6) 



dt 



Differentiating (|5.6|) gives 



d_( d_\ _ (Pu d 2 u d i 

dt\dt U ~ dt* + dxWtdt X [ } 



= w—a + a—w. 
dt dt 



Now, we can either compute directly or use equation (|2.8p . S = 0, a = eH to see 
that 

Jt U = ° H - (5 ' 8) 
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In turn, (15. 7| becomes 
dt 2 



eHw + a—w + awi 
at 



d 



d 2 u \ d 



dx l dt dt 



eHw + cr—w H 



dt ■ ' w \d^dt 



Since w, = i 8 kl UkUu and Jjtu = ^ ^ u jS^m.i we obtain 



w 3 dx l dt ' 



9i 2 w dx l dt 



3 w 2 



S M S ij 



UklUiU-j. 



In addition, from (15. we have 



and from 



1 

G= 2 



to 1 + LDit 



1, 2 x 1 



5^ + »)-5 T 



(5.9) 



(5.10) 



Moreover, the second fundamental form is given by 



and the induced metric and its inverse by 



9ij = % + UiUj, g lJ = 8 %3 - -\s ik 6 jl u k ui, 

w 2 



so that 



1 



Hw = g l3 u %3 = IS* 3 -S'H^UkUi u 



(5.11) 



This leads us to the equation (|5.1[) for graphs over a flat subspacc Z 1 - . 

Step 2: Hyperbolicity of the second-order equation. The equation (|5.1|) is hyperbolic 

if and only if the following matrix 



( -1 2-U! 

f ux 

in ± 



a,f3=0,...,n 



(A* 3 ). . , 



with 



satisfies 



A« := (e + cr 2 ).g y ' - cr 2 ^' , », j = 1, . . . , n 



d k :=det(( J 4^) ( 



< 0, < k < 



I ot,P=0,...,kj 

Namely, fix any point p £ Z , and choose an orthonormal basis e\, . . . , e n spanning 
Z 1 - such that at p 



Du(p) — u\e\, ui — 0, 2 < i < n. 
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In this basis the matrix A at p takes the (symmetric) form 



( Aa ) a ,B=0,. 



■1 



- ui 








V 







o\ 






e 



This gives in view of u 2 = |_Du| 2 

d = -1, d k 



< 



k> 1. 



Since p was arbitrary, this proves that equation (|5.ip is indeed hyperbolic. 
Step 3: First-order formulation. By introducing the first-order variables a = u t /w 
and bi :— Ui, and regarding w as a function of b = (bi), we can rewrite (I5.ip in the 
form 



da bi 3cr 56, 
Ot-^S'^-e(a) 

at 



w(b) dxi 



w(b) dxi 
w(b) t^— - a—?^-$ -z—r = 0. 



0, 



dx r w(b) dx 1 

This is a first order nonlinear system in <r, b, which can be checked to take the 
desired conservative form (|5.2[) . 

Considering the function E = (<r, b) introduced in the theorem, we also compute 

8 2 E 



8a 2 
8 2 E 

dadbk 
d 2 E 

db.dbk 



w, 



= w 1 abi. 



w 



-.I ( ^ + „,( { ,-^) 



which is a non-negative matrix, since for all scalar Y and vector X = (X^) (not 
both zero) the Hessian evaluated at (Y, X) equals 



X^bjX k b k 
w 2 



or equivalently 



w(Y + v^X j Y + w ~ 3 l° 2 {{X j X k 5 jk ) - 2{b l X°) 2 ^j 
+ vT^XixHjk + ™- 3 ^(a 2 + n) ((X^ X h 8 jk ){b i b l 8 u ) - (X*bj) a ) . 

This expression is positive if and only if we impose the restriction bibj5 tJ < 1/2. 
(Each term in the above decomposition has a positive sign.) 
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Since the equation (15. 2p has the form of a system of conservation laws and admits 
a convex entropy, it can be put in a symmetric hyperbolic form. Indeed, introducing 
the variables 

a := u t = <rw(b), Cj := e(er)- 



w(b) : 



0. 



(5.12) 



which is nothing but the gradient of E, we obtain 

Then, by expressing (implicitly) a and &j as functions a, bi of the new unknowns a 
and Ci, one can check that the above system is symmetric, in the sense that 

do db~i . 

_(«,<.) = _(«,<.). 

In turn, the system is locally well-posed in H s with s > 1 + n/2. □ 

Weak solutions to the normal flow of graphs. To define weak solutions we 
rely on the conservative form exhibited in (|5.2[) . 

Definition 5.2. A Lipschitz continuous map u : [0,T) x 1" — > R is called a weak 
solution to the HMCF' equation for graphs (|5.1|) if and only if for every test-function 
9 : [0, T) x R" -> R (that is, compactly supported C°° functions) 

2 



(0.T) 

where w(Du) 2 
inequality 



du 89 
~dtdt 



fiij du du 




(w(Du))~ 



du 
~dt 



du 89 



n 



S ij (w(Du))" 1 dxdt = 0. 



dx l dxi 

Tjx I ~dxir ^ ^ s ca H e d a n entropy solution if, moreover, the 



(w(Du)Y 



du 
di 



n w{Du)9dxdt < 



for every non-negative test-function 9. 

We have the following uniqueness result, which relies on the fact that the energy 
is convex in the variables tr, b. 

Theorem 5.3 (Uniqueness of classical solutions within the class of entropy solu- 
tions). Given e > 0, there exists a constant C e such that the following property 
holds. Let u be a Lipschitz continuous entropy solution and u' be a solution of 
class C 2 , both being defined up to some time T > and satisfying the uniform 
hyperbolicity condition 



\Du\ z < 



1 



\Du'\ 2 < 



1 - e 



2 ' 2 
Then, provided the Lipschitz norm of u if less than C £ and the C 2 norm of u' is 
less than C' e then for all times t € [0, T) 

du 
~di 



< C t e 



c t t 



du' 


2 






+ 


\Du- 


" ~dt 






du 


du' 


/..( 


di 


~ ~dt 



(t, x) dx 



\Du-Du'\ (t,x)dx. 
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Proof. Under the assumptions made in the theorem, consider the expression 
Q = Q(t, x) = E(a, b) - E{a'.b') - j^^(a', b')((a, u) - (a', u')) 

= ((1 + \Du\ 2 )- Vt| 2 + nj Vl + \Du\ 2 - ((1 + iDu'f)- 1 ]^ 2 + n) V 1 + \ D ^'? 

and note that for some constant C\ > 

(|n t - u' t \ 2 + \Du - Du'\ 2 ) <Q<C\ flu, - u' t \ 2 + \Du - Du'\ 2 ). 

On the other hand, a direct calculation using the fact that u is an entropy solution 
and u' is a classical solution yields the inequality 

d f f 

1 Q(t,x)dx <C 2 \ Q{t,x)dx, 



dt 

where the constant C 2 depends upon up to second-order derivatives of the solu- 
tion v! . The conclusion follows from Gronwall's inequality and the fact that Q is 
comparable with | ut — u' t \ 2 + \ Du — Dv! | 2 . □ 



Global existence of one-dimensional graphs. In view of Theorem 15.11 the 
hyperbolic mean curvature flow equation in the case n = 1 reads 

If U 2 \ U xx -UxUtUxt U x U 2 t U xx 

u tt = 77 7— o + 1 7— o + 2- 



u tt = : „. 2 , 2 u xx + 2 2 u^. (5.13) 



2V1 + U 2 J l + ul l + U 2 x (l + lt 2 ,) 2 

= 7777 77^" ( u t + 1 + u x ~ 2u_.lt, ) + 2- — . 

2(1 + 1*2)2 \ t x x t) l + u\ 

Therefore, we find 

1 + u\ + u 2 - 2u 2 x u 2 u x u t 

2(1+ U 2)2 UXX + Z 1 +U 

This equation written on the real line with initial data 

u(x, 0) = uo(x), ut(0,x) = u%(x), 

describes the vibrations of an infinitely long string, with initial position uq and 
initial velocity u\. 

Relying on the definition of weak solutions introduced earlier for general dimen- 
sions, we now prove: 

Theorem 5.4 (Global existence of weak solutions). There exists a constant Sq > 
such that given any initial data Uq,ui : R — > M such that 

TV(u Q . x ) + TV( Ul ) <5 Q , 

the initial-value problem for the equation (HMC'F') in the second- order form (|5.13[) 
admits an entropy solution u — u(t, x) such that the functions u x and u x have 
bounded variation in space, uniformly in time. 

Proof. Since 

C~ 1 \ 1 -\- u 2 -\r u 2 e 

u x u t l+ul+u 2 t -2ulu 2 t I = 2(1 + U 2 ) 2 = ~~ V? < °' 

this equation is always hyperbolic. We introduce the variables 

a := — , o := u x . 

w 
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These are both conservative quantities. For b this is trivial since 



h = u xt = {u t ) x = (aVT + b 2 ) x . (5.14) 

d_ 

c)u 



For a we may either see this directly from the conservation law (|4.2| with Y — 



or we may compute 

u tt u x u t u xt 



w 



l + u 2 x + u 2 - 2u 2 x u\ n u x u t u x u t u xt 
u xx + 2 — — u x t 3 — 

u x (u?) 





2w 5 




1 - 


hul + u 2 - 


2u 2 x u 2 




2w 5 




1 - 


i-ul + u 2 - 


2u 2 x u 2 




2w 5 





2w 3 

2 ( u x \ 



thus 



a, 



2w 3 



2w 3 V 2w3 

u x ^u x u 2 \ ( u x (\ + u 2 + u 2 )\ ( (l + a 2 ) b 



2 

^xx { Hx^t 



2w 2w 3 J x V 2 ^ 3 Jx \2yTTF 
This can be rewritten in the form 

ab 1 + a 2 

at = vTW ax + 2(1 + 62)3/2 b * 



and equation (|5.14[) is equivalent to 



ab 



b t = \/l + b 2 a x + - b x . 

Vl + b 2 

Combining the last two equations gives the system 

Q -A(a,b) Q =0, (5.15) 

where 

( b 1 + fl2 \ 

A(a,6) = -J= 2(1 + 6 2 ) . (5.16) 

^ /r +^l v l + 6 2 ab J 

Recall that a necessary and sufficient condition for a quantity r\ — r)(a,b) to be 
a conserved quantity is for 



T]aa 'Hab 



A 



y Vab Vbb / 

to be a symmetric matrix. For the nonlinear hyperbolic system under consideration 
this gives 

f-, .2s 1 + a 2 , 

abr/ab + (1 + b ) r] bi = — ^ i] aa + abr] ab , 

that is, 

2(! + fe2 ) 2 

Vaa = — 7—, 5— Vbb- (5-17) 

1 + ar 
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This is clearly a linear hyperbolic equation. From this we compute for example 
that ab — ^^j^ 1 is a conserved quantity, a fact which also follows directly from (14. 2|) 
with Y = 

ox 

As one easily computes, the two eigenvalues of A are given by 



VTT¥ 

while the eigenspaces are spanned by the vectors 



A± = ^=^ U6±J±±M, (5.18) 



± 



v /5 a+ b2 ) j . (5.19) 



Regarded as a function of a and b the gradient of A± equals 

/ b± 



DX± = ' 




(5.20) 



hence 



I / ± VT+q^ / g \ a 6 /l + a 2 



^/TTF^ V2(i + fo 2 ) v 1 + 62 1 + 62 

3a 3ut 

= 2(1 + fo2)3/2 = 2^4- 

Hence, the hyperbolic system under consideration is not genuinely nonlinear in the 
sense of Lax. 

However, we observe that the genuine nonlinearity is lost on a hypersurface (that 
is, {ut — 0}) which is itself non-degenerate, in the sense that 

(D(DX± , /i±), (i±) along the hypersurface u t = 0. 

Therefore, we are in a position to apply to the system of conservation laws 

I 2VTTF ) 



{ay/l + b 2 ) x = 0, 



the global existence theorem in Iguchi and LeFloch [4] , which provides the existence 
of a solution with bounded variation when the initial data have small bounded 
variation. □ 



6. Local-in-time existence 

We now turn to the discussion of the existence of solutions to the normal hy- 
perbolic mean curvature flow (HMCF'), where M t cannot necessarily be written as 
an entire graph over a flat subspace. We use standard notation and, in particular, 
denote by H£ oc (M) the Sobolev space of locally squared integrable (tensor- valued) 
maps defined on M whose all s-order derivatives (in one local chart and in the 
distributional sense) are also locally squared integrable. 
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Theorem 6.1. Let M be a smooth, orientable compact manifold with dimension 
n, and F : M — > R n+1 be an immersion of M in the Euclidian space. Given 
a (scalar) normal velocity field 7T : M — > M in the Sobolev space H s+l {M) with 
s > 1 + n/2, there exists a unique flow F : [0,T) x M — > R n+1 in the space 
L°°([0,T), H s+1 (M))n Lip([0,T), H S (M)) which is defined on some maximal time 
interval and satisfies the normal hyperbolic mean curvature flow equation (HMCF'), 
together with the initial conditions 

F(0) = F, ^(0) = <f 1/(0). 

We will provide two different arguments to handle the equation (HMCF'). 

Let us us cover the manifold with finitely many local charts, chosen in such a 
way that the manifold can be viewed locally as a graph over its tangent plane at 
some point. In each local chart, we apply the local existence theorem for graphs 
established in Thcorem l5.ll Indeed, due to the property of finite speed of propaga- 
tion satisfied by hyperbolic equations, all of the arguments therein can be localized 
in space and apply in each coordinate patch. Then, by patching together these 
local solutions and using the fact that only finitely many charts suffice to cover 
the manifold M, we can find a sufficiently small T such that every local solution 
is defined within the time interval [0, T), at least. This completes the proof of the 
theorem. 

The rest of this section is devoted to provide a second proof of Theorem l6 . 1 1 which 
is also of interest in its own sake. We will now express (HMCF') as a single scalar 
equation in terms of a height function u with respect to a fixed initial hypersurface. 
To this end let us discuss the case of flows F : [0, T) x M — > M" +1 such that each 
M t := F(t, M) is an entire graph over a fixed reference manifold £ given by an 
immersion G : M — > S C M™ +1 . If each M t can be written as a graph over E, there 
must exist a family of smooth height functions u : [0, T) X M — ► M and a family 
£ : [0, T) X M — > M of diffeomorphisms such that 

F(t,x) - G(£(t,x)) +u(t,Z(t,x))n(Z(t,x)), (6.1) 

where n is the inward unit normal along E. Let us denote the metric on M = £ 
by <Jij and the second fundamental form by The induced connection on M with 
respect to a will be denoted by D. 
The tangent vectors take the form 

Fi(t,x) - (G i (C(t,x))+u i (t,C(*,x))n(C(t,a;)) 

-u(t, £(t, x))t* m x))G k (£(*, x)))£ (t, x) 

where in this section a raised index will be raised with respect to the metric er, i.e. 
T k _ (jfe'Ty, it w iH be convenient to define the following tensor 

Nij := Gij - unj. 

Then the tangent vectors can be written in the form 

Fi = (ujH + NjGi) £j (6.2) 

and for the second derivative i 7 ^ := d ®ig x j we get 
F 3 = {{u kl + Nl n r rnl ) n + (DiN™ - u k r?) G m ) g$ + {u k n + N™G m ) 
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The induced metric tensor gtj = (Fi, Fj) is 



<to- = (ufcUi+iWVH£(£$- (6-3) 

In the following we will assume that u is sufficiently small, so that the symmetric 
tensor is invertible and we denote its inverse by . Let us define 

to := \Jl + N k N a u k ui. 
The inward unit normal along M t is then determined by 



v= - (n~N kl u k Gi) , (6.4) 

so that 

hij = (Fij,i/) 



i (u kl + N^Tml + N rm U r (u k T ml - A^fcrn)) (6.5) 



We need expressions for -^F and 4&F. From (|6.ip we obtain 



where the subscript t in u t denotes a partial derivative with respect to t, i.e. Ut 

du 
dt 



~ . Rearranging terms gives 



iF=(u t +u k ^-)n + Ni^ Gl . (6.6) 



dt V dt J dt 

This implies the relations 



and 



dt w 



s '' = (i F '^) = (^+^^)-' ,(is » 



where 

9ki ■= u k u t + N^N ml . 
We differentiate (|6.6|l with respect to time and compute 

d? F = { Utt + 2utk ^r +Uk ^) n -{ Ut+Uk ^r) T ^ Gl 

Hence, we have 

~dt? = { U - + 2U ^ + U ^ + N ^ll) n (6 ' 9) 
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and, consequently, 



1 / d£ k d 2 £ k d£ k dP 



Therefore, we have 



1 / d£ k d£ k df* 



- AWi iv r „„f f - + L + »| n'iVr« m f (6.10) 



and, moreover, 



( dt k d 2 £ k d£ k d£ s \ 



Reordering gives the final formula 



d$, k d 2 i k . , ^d£ s df fc 

UjUu + 2ujU tk — + 9jk-j^- + (UjN k T sl + N ol D s N k)-^-i^ 

^nr l s[u t + u k d p i d p i ii. (6-1 1 ) 



dt J dt 

We summarize our results in the following proposition. 

Proposition 6.2. If F : [0, T) x M — > M" +1 is an arbitrary flow, where each M t — 
F(t,M) is represented as a graph over G : M — > S as above, then (u, (£, k )k=i,...,n) 
is a solution of the coupled system 

a = (6.12) 
£ = (6-13) 



a 



1 / d£ fc (if fc d£* 



</c </c , ( ae\ # 



" A «"^1r + h + ) (6 ' 14) 



(df k d 2 £ k d£ s d£ k 

Uj u tt + 2u jUtk ^- + ~9jk-^r + i u i N kTsl + N H D ' N k)-£-£ 



-N^[u t + u k ^ y%)$. (6.15) 
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Note, that up to this point we have not chosen a particular flow. This will be 
done in the next step. For (HMCF'), we have Si = 0, a = eH, A = —Vie. From 
(|6.13[) we obtain 

-|- = -g jk u jUt) (6.16) 
where g 3 denotes the inverse of gjk- Inserting this into (|6.14|) we get 

cH = - {u tt - 2g jk u 3 u t utk + L] , (6.17) 

where L is a term of first order in u. This, by the compactness of £, clearly is 
a uniformly hyperbolic equation, provided u is sufficiently small in C 1 (E). So we 
can provide a short-time solution for (HMCF'), if we assume uu—q = 0, i.e. if 

£ = Mo. After having solved (|6.17j) we can solve (|6.16|) for ^fr-. If u is sufficiently 

small in C 1 (E), then ^4r- will be small as well so that by choosing £ t= o = Ids we 
obtain a family of diffeomorphisms £(t) solving (|6.15[) . This completes the proof of 
Theorem 16.11 

7. Finite-time blow-up results 

An example of finite-time blow-up. The hyperbolic mean curvature flow may 
blow-up in finite time in a way that is completely analogous to the standard mean 
curvature flow. We provide here a typical example. 

Let F Q : S n -> R n+1 be a round sphere of radius r . If F : [0, T) x S n R" +1 
is a solution of (HMCF') with initial data F(0, x) = F (x) and ^F(0, x) = <tqu(x) 
with a constant <7o G R> then M t := F(t, S n ) is a concentric sphere with radius r(t). 
In this case, the hyperbolic mean curvature flow reduces to the ordinary differential 
equation (ODE) 

n, n 2 
rr + -(rf + — = 0, 

r(0) = r , 
r(0) = -* Q - 

This second-order equation can be reduced to the following first-order ODE 



T (7.1) 
- J(n + ctq) (-7) -n, if cr > 0, 



r(0) = ro. 



The solution depends upon the dimension n. In the case n = 1, we obtain the 
cycloid 



c / c 

1 + c arctan t 1 

r V r 



—t — 7*0(7 — c arctan (TO) if 00 < 0, 
t + 7'o(To + c arctan (Tq, if (To > 0, 
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where c = ro(l + <7q). 

On the other hand, in the case n = 2, we obtain the explicit solution 



r(t) 



•2 - 2r a t - 2t 2 . 



If (To < 0, the sphere begins to expand until it starts to shrink and eventually 
collapses to a point in a finite time T, given by 



T = y { - ao 



ol + 2 ) . 



In the above situation, one can avoid the formation of singularities by rescaling 
the metric according to its volume. 

Blow-up estimates based on the mean and total mean curvature. In some 
situations it is possible to derive blow-up results from the behavior of the mean or 
total mean curvature of the system. To this end let us define the function 

7(c) := -^= arctan ( -^= ) , 
Vn \y/nj 



which satisfies 



and 



7' = -- 1 



n 21 + 1 



d ,d 

dt 1 = 1 dt a = H - 

Given U C M we define 

fu{t) ■= j(a + 7 e)dn u E v := J ed\x t - 
u u 

Note that -^(edfit) = implies, that Ejj does not depend on t. 

Proposition 7.1. For any integrable U C M and any < t\ < t 2 < T one has 



Hd^itdt 



ti u 



and 



J J Hcd[i t dt 

ti u 



< 



^fuM-fvit^^Eu 



fu(t 2 ) - fu(h) + J o-dfj, tl - J o-d^ t2 

u u 

2(tt+1) 



(7.2) 



Ejj . 



(7.3) 



In particular, if T = oo, then for any e > and any choice of integrable U C M 
there exists a sequence tk — > oo such that 



Hdfi tk 



u 



< e, 



fc e N 



THE HYPERBOLIC MEAN CURVATURE FLOW 



25 



Hedn tk 



< e, keN. 



Proof. Since J^cx = eH and 4r d/it — —crHdfit we compute 

^fu(t) = J ( cH + Hc + l acH ~ ( a + 7e)o-H)dfM = n J Hdfjt t . 
u u 
The function — + 7 is a monotone increasing function in a and 



7T G 7T 
7= < - +7 < ~F=- 

Jn e Jn 



Therefore the function 
satisfies 

Then we obtain 



fu{t) 



- + 7 edfi t 
e 



v 



^=E u <f u {t) <^=Eu. 
yjn Jn 



(7.4) 



dt 



fu(t) =n Hdn t 



which implies 1|7.2[) . From the observation 

< J (a± Jn) 2 du t = J{cr 2 +n± 2Jna)du t = 2E V ± 2Vn J adu t 



we conclude that 



adut 



< 



Eu 



(7.5) 



Moreover, we have 
d 
dt 

and thus 



J a dfit = j "( e ~ cr 2 )Hdfj, t = — J cHdut + n J Hdut 



j{fu{t)-Jad^ 
\ u 



= I eHdfif 



This and (EH), (USD imply fTS]) . 

We can do even better: 
Proposition 7.2. For any x £ M and any < t\ < t% < T one has 

«2 



□ 



H(x,t)dt 



2tt 



In particular, if T = 00. then for any e > and any x € M there exists a sequence 
tk — » 00 such that 

\H(x,t k )\<e, fceN. 
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Proof. This follows directly by integrating 
and from 

M < 4=- (7.6) 

y/n 

a 

Proposition 7.3. Consider the flow associated with a closed curve Cc I 2 with 
non-vanishing rotation number x(C). Then T < oo. 

Proof. The rotation number of a curve is given by 



Hdii. 



This is a topological invariant, hence in particular 4iX(pt) — f° r an smooth 
deformations Ct of C. It follows that 

Hdfj, t 

is a (non-zero) constant. Hence, by Proposition 17. II we must have T < oo. □ 
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